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It has been pointed out that the recent BaBar data on the ivy* — > 7 transition form factor 
F 1T1 {Q 2 ) at low (high) momentum transfer squared Q 2 indicate an asymptotic (flat) pion distribution 
I amplitude. These seemingly contradictory observations can be reconciled in the kr factorization 

^vj theorem: the increase of the measured Q 2 F 1T1 {Q 2 ) for Q 2 > 10 GeV 2 is explained by convoluting 

a kr dependent hard kernel with a flat pion distribution amplitude, kx being a parton transverse 
momentum. The low Q 2 data are accommodated by including the resummation of a s In 2 x, x being 
a parton momentum fraction, which provides a stronger suppression at the endpoints of x. The 
next-to-leading-order correction to the pion transition form factor is found to be less than 20% in 
the considered range of Q 2 . 

' <-]', PACS numbers: 12.38.Bx, 12.39.St, 13.40.Gp 

Oh| The asymptotic and soft behaviors of the pion transition form factor F V1 {Q 2 ) involved in the process 7T7* — * 7 

have been derived [l|: 

lim Q 2 F^{Q 2 ) = V2U = 0.185 , (1) 

Q 2 ^oc 

: J* F^(Q 2 ) = -^L , (2) 
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Q 2 being the momentum transfer squared carried by the virtual photon, and f„ = 0.131 GeV the pion decay 
constant. The former is predicted by perturbative QCD (PQCD) in the collinear factorization theorem 0, 0, 0, 
HI HQ' while the latter is determined from the axial anomaly in the chiral limit. However, the recent BaBar 
data on F 7r7 (Q 2 ) exhibits an intriguing dependence on Q 2 Q 2 F 7ri (Q 2 ) exceeds Eq. ([]]) for Q 2 > 10 GeV 2 , 
and continues to grow up to Q 2 « 40 GeV 2 as shown in Fig. [TJ It has been commented 0] that this behavior 
cannot be explained by perturbative effects, like higher-order or higher-twist contributions. 

It is known that one can extract nonperturbative information on the shape of the leading-twist pion distribu- 
tion amplitude (DA) from the measurement of F 1ll {Q 2 ) Q. Inspired by the BaBar data, a model for the pion 
' DA, differing from those investigated in R|, has been proposed fioj . 

<j> n (x, t i ) = N+(l-N)6x(l-x), (3) 

with N being a free constant. This model remains finite at the endpoints of the momentum fraction x — > 0, 1 
for the normalization point un = 0.6 ~ 0.8 GeV. It was argued [lfj that the pion DAs from the instanton 
theory of QCD vacuum [HI, [ij], from the the Nambu-J ona-Lasinio model |13[, from the chiral quark model 
[l4| . and from the large- N c Regge model il5j are expected to be rather fiat. At the same time, the hard kernel, 
proportional to the internal quark propagator, was modified by introducing an infrared regulator m 2 [Tol |. Note 
that an endpoint singularity would be developed for the model in Eq. ([3]) without the infrared regulator. Then 
the factorization formula 

3 J x + m z /Q' 1 

leading to a logarithmic increase InQ 2 , explains the growth of Q 2 F ni (Q 2 ) observed by BaBar with the tuned 
parameters N — 1.3 ± 0.2 and m — 0.65 ± 0.05 GeV. The regulator m 2 was interpreted as the inverse instanton 
size, which sets the scale for nonperturbative effects in a quark propagator [l(| ■ 
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A similar conclusion on the shape of the leading-twist pion DA has been drawn from the BaBar data in [16(] , 
but with the difference in the formalism that the parton transverse momentum kx was taken into account. The 
hard kernel is the same as in the collinear factorization, and the In Q 2 dependence comes from the integration 
of the pion wave function over the intrinsic parton kx up to xQ [IB]: 



Q 2 F V ~ I {Q 2 ) 



V2f„ f 1 dx 



Gtt 



xQ 



kxdkx ipTi(x, kx) 



(5) 



with the model 



ip„(x, k T ) = 2 ^^( x ) exp 
x(l — x)a 



2ax(l — x) 



(G) 



The cutoff of kx at xQ guarantees that the integration over x down to is finite even for a flat distribution 
4>tt(x), with which the BaBar data are accommodated by choosing the parameter a = 0.53 GeV 2 . A concern 
was raised for a broad pion DA in the conventional collinear factorization [1(1 : the next-to- leading-order (NLO) 
correction to the hard kernel [l7j is huge for the renormalization scale /i ~ Q, or /i has to be much smaller than 
the scale Aqcd in order to diminish the NLO correction. For the calculation of the next-to-next-to-leading-order 
correction to the pion transition form factor in the collinear factorization, refer to [HI]. 

The BaBar data are consistent with those of CELLO [l9| and CLEO 12011 below Q 2 w 9 GeV 2 , from which 
an endpoint suppressed shape of the pion DA has been extracted [2 j| . [2 2l|. Taking into account high-power 
corrections, the asymptotic pion DA can also explain the CLEO data [23|. The QCD sum rule analyses of the 
pion DA and of the pion transition form factor led to an endpoint suppressed model 0, Hi], [U, [2^, [26|, [27], [28[ . 
It has been noticed [2!| that the CLEO data favor the asymptotic form rather than the Chernyak-Zhitnitsky 
one Q, which emphasizes the endpoint region. The E791 di-jets data support the nearly asymptotic pion DA 
above M 2 — 10 GeV 2 [30], Mj being the mass squared of the di-jets. A review on the determination of the 
leading- twist pion DA can be found in (3l| . It seems that the measurements of the pion transition form factor 
at low and high Q 2 imply different shapes of the pion DA, which are contradictory to each other. 

In this letter we shall study the process ttj* — > 7 using the kx factorization theorem [H, [H, [H, HH, [H, H3] ■ 
If a flat pion DA is favored, it will enhance the contribution from the region with small momentum fraction x. 
Once the small x region is important, the parton transverse momentum squared k T , being of the same order as 
xQ 2 in the internal quark propagator, is not negligible. The regulator m 2 in Eq. ^ can be interpreted as k T 
at LO, whose inclusion smears the endpoint singularity from x — > 0. Besides, the loop correction to the virtual 
photon vertex generates the large double logarithms a s ln 2 (Q 2 /fc^) and a s In 2 x in the small x region [3a ]. The 
former is absorbed into the kx dependent wave function, and organized to all orders by the kx resummation 
35 , 38 , 40( . The latter is absorbed into a jet function, and organized to all orders by the threshold resummation 
4U l42l |43| ] . With the above resummation factors, a flat pion DA does not develop an endpoint singularity. 

The NLO correction to the hard kernel of the pion transition form factor has been computed in the kx 
factorization [38[. Fourier transforming Eq. (44) in [3a] into the impact parameter bx space, we derive 



F^{Q 2 ) = [ dx [ b T db T Mx)eM-S{x,bT,Q,^}S t (x,Q)K (^Qbx) 

Jo Jo 



1 - 



4tt 



-C f In 



_ , 7r + 2 In a; + 3 



(7) 



where Cf = 4/3 is a color factor and 7^ the Euler constant. The NLO term in the square brackets implies 
that the renormalization scale /i should be chosen as fj? ~ Oi^fxQ /bx) in order to minimize the logarithm. 
However, we have always set /i = max(y/xQ , 1/bx) in our previous analysis of exclusive processes [3a], so we 
shall continue to adopt this choice here. It has been confirmed that the above two choices of \x produce almost 
identical results. The Sudakov exponent from the kx resummation is given by 



S(x, b T , Q,n) = s f x ^=> b T 



{l-x)-j=,bx 



^7 ? (a s (/x)), 

l/b T A* 



(8) 



where the explicit expression of the function s and the quark anomalous dimension j q can be found in [44l . |4£ 
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g (GeV) 2 



1.7 0.75 0.47 0.32 0.25 



TABLE I: The power c in the threshold resummation at different Q 



The power c « 0.3 in the parametrization St(x,Q) for the threshold resummation [46 



(9) 



was derived at the scale of the B meson mass ms = 5.28 GeV. The power-law behavior should be modified 
as Q 2 runs within a large range in the present case. We determine the Q 2 dependence of the parameter c by 
repeating the procedure in Appendix D of [4^ | , which involves the best fit of Eq. (|9|) to the exact resummation 
formula in the Mellin space N 



S t (N,Q) = exp 



- l ' N dz 



with the anomalous dimension 



IK = C F + 

7T \ 7T 



(1-z) 



VAQ72 



67 ^ 2 \ 


5 


36 12 J 


18 



(10) 



(11) 



rif being the number of quark flavors and Ca = 3 a color factor. The outcome is listed in Table [l] where the 
rapid increase of c with decreasing Q arises from the exponentiated radiative correction proportional to a s (Q 2 ). 
It becomes more difficult to determine c for Q < 3 GeV, because the allowed interval of A" shrinks. Therefore, 
we freeze c at c = 1, when it exceeds unity. To simplify the analysis, we propose a parabolic parametrization 



c -> c(Q 2 ) = 0.04Q 2 - 0.51Q + 1.87 



(12) 



if c < 1. A large value c ~ 1, resulting in a quick falloff as x — > 0, will improve the agreement of our prediction 
with the low Q 2 data, when the flat model is adopted. That is, a careful treatment of the threshold resummation 
effect with a running c(Q 2 ) is crucial for accommodating the BaBar data in both the low and high Q 2 regions. 
We adopt Eq. ([6]) for the pion wave function, which is written as 



(f>Tt{x,b T ) = <Mx)exp 



-ijx(1 



(13) 



in the impact parameter space. The Gaussian form in b 2 -. is the same as proposed in |29j . and consistent with 
that from [47j. For the asymptotic model, ^{x) — 6x(l — x), the parameter a can be fixed by the normalization 
condition 48] 



dx I d b T (j} v {x,b T ) = —p: , 



(14) 



giving a = 4ir 2 f 2 = 0.677 GeV 2 . For the flat model, (p^(x) — 1, the above normalization condition does not 
apply, and the parameter a will be tuned to fit the BaBar data. In the analysis below we choose a — 2.5 GeV 2 . 
Note that the renormalization-group evolution of the Gegenbauer expansion starting with the asymptotic form 
has been known. However, the evolution starting with the flat form has not yet been studied. Hence, the models 
adopted here should be regarded as being defined at a low normalization point /x . 

The LO and NLO predictions for the pion transition form factor from the fcy factorization are presented in 
Fig. [U It is observed that the curves for Q 2 F^ 1 (Q 2 ) from the asymptotic pion DA saturates quickly for Q 2 > 5 
GeV 2 . That is, it is difficult to explain the continuous growth of Q 2 F 7ri (Q 2 ) above Q 2 ps 10 GeV 2 using the 
asymptotic model. The overshoot in the low Q 2 region is expected, since this portion of data favors an endpoint 
suppressed shape [2ll I22I ]. Note that the overshoot can be resolved within the theoretical uncertainty from 
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FIG. 1: Q 2 dependence of Q 2 F nj (Q 2 ) up to LO and NLO from the kr factorization for the asymptotic (red, dotted and 
dot-dashed lines, respectively) and flat (blue, dashed and solid lines, respectively) pion DAs. The points with errors are 
the data from CELLO pj], CLEO [13], and BaBar [|. 




CT [GeV^] 

FIG. 2: Q 2 dependence of Q 2 Ft I1 {Q 2 ) up to NLO with the power c = 0.3 (red, dashed line) and c(Q 2 ) (blue, solid line) 
in Eq. dJ). 



varying parameters in our formalism. For example, multiplying the arguments xQ/y/2 and (1 — x)Q/y/2 by a 
factor 2 in Eq. ((8]) [35| leads to stronger suppression at x — > 0, and better agreement with the data at low Q 2 . 
The curves from the flat pion DA show a good fit to the BaBar data, whose logarithmic increase with Q 2 is a 
combined consequence of the inclusion of the parton kr and the employment of the flat model. We stress that 
the fit would deteriorate in the low Q 2 region without a careful treatment of the threshold resummation effect. 
To highlight this effect, we compare the results from the flat model with a constant power c = 0.3 and with a 
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FIG. 3: Q 2 dependence of Q 2 F 7rl {Q 2 ) from (a) the fiat model with different a. and (b) Eq. (|I5[) with different a-i- 



running c(Q 2 ) in Fig. O It is found that the former exceeds the data around Q 2 w 10 GeV 2 . The results from 
the asymptotic model is almost independent of the running of c(Q 2 ), because the asymptotic model decreases 
quickly enough at small x by itself. The NLO correction, amounting up to 20% for the flat pion DA as shown 
in Fig. [1] is under control in the kx factorization. This observation is opposite to that postulated in [lOj , and 
to that obtained in the collinear factorization [iH ]. 

We should keep in mind that the BaBar data suffer from a large uncertainty in the high Q 2 region. In 
particular, the data point for Q 2 = 27.31 GeV 2 drops, and is consistent with the asymptotic limit in Eq. (JTJ) . 
Hence, we examine the dependence of our prediction on the parameter a associated with the intrinsic bx 
dependence in Eq. (|13p . The results are displayed in Fig. [3)Ja), indicating that the slope of the curves is 
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insensitive to the variation of a especially in the high Q 2 region. If including the higher Gegenbauer term into 
the asymptotic model 

4>v{x) = 6x(l - x) [l + a 2 C 2 /2 (2x - 1)1 , (15) 

with the Gegenbauer polynomial C 2 ^ 2 (t) = (3/2)(5i 2 — 1), the slope of the curves changes significantly with a 2 
as shown in Fig. [3jb). The prediction is enhanced by 10% for a 2 = 0.2 in the high Q 2 region, and saturates 
around Q 2 10 GeV 2 . Only with an extremely large a 2 ~ 0.8, which also emphasizes the endpoint regions, 
the prediction becomes close to that from the flat pion DA. It has been mentioned [10j that a very broad pion 
DA, vanishing at the endpoints, but with a rapid increase (</4(0)/6 S> 1), serves the purpose of accommodating 
the BaBar data. However, the big value a 2 « 0.8 may render the Gegenbauer expansion starting with the 
asymptotic model questionable, as claimed in [13] • In summary, the precise measurement of the pion transition 
form factor in the high Q 2 region will settle down the issue on the shape of the leading-twist pion DA: if the 
growth with Q 2 observed currently persists, the flat model is favored. If the growth becomes milder in the 
future, the asymptotic model with nonvanishing higher Gegenbauer terms will be still allowed. 

Finally, we comment on the kx factorization formalism developed in [42| , which involves on-shell partons in 
the calculation of hard kernels. The usual procedure of factorizing a kx dependent wave function is to neglect 
the minus component k~ in a hard kernel, and then to integrate out k~ in the general wave function that 
depends on the four components k + = xp + , k~~ , and kx, 

tp(x,k T ) = I dk~^(x,k~,kr). (16) 

This is the reason why a parton, participating in the hard scattering, carries an off-shell momentum (fc + , 0, kx) 
in our formalism 49]. The prescription for calculating a kx dependent hard kernel, the definition for a kx 
dependent wave function, and the gauge invariance of the kx factorization have been carefully addressed in 
[38L [49| , which confront the criticism from [50L l5lT ] 1 . The transverse momentum squared k T appears in the hard 
kernel for the pion transition form factor through the internal quark invariant mass, (k — -P 7 ) 2 = —xQ 2 — k T , 
for the final-state photon momentum P 1 = (0, P~,0x) and k 2 = —k T . Assuming on-shell partons, the hard 
kernel is independent of k T because of (k — P 7 ) 2 = —xQ 2 for k 2 — 0. Therefore, it will be difficult to explain 
the BaBar data using the formalism of [421 ]. 

We thank Ahmed Ali for useful comments. This work was supported by the National Center for Theoretical 
Sciences and National Science Council of R.O.C. under Grant No. NSC-95-2112-M-050-MY3. 
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